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Abstract. —Phylogenetic relationships among 20 currently recognized species of Botrychium sub¬ 
genus Botrychium were investigated using nucleotide sequences from the plastid gene rbcL. Anal¬ 
ysis of eight diploid species produced a single most parsimonious tree with three distinct c .ades. 
The lanceolatum clade (one species) was sister to all other species. The lunaria clade (two species) 
was sister to the simplex-campestre clade (five species) which contained two weakly supported 
subclades. Single populations of B. lunaria and B. simplex were more closely related to other 
species in the same geographic region than to populations of their own species from distant geo¬ 
graphic regions. In a second analysis that included both diploid and polyploid species, eight of 
twelve polyploid taxa had rbcL sequences identical to those of known diploids, which suggested 
relatively recent origins for these eight polyploid species. Four polyploids, B. ascendens, B. par- 
adoxum, B. Xwatertonense, and B. minganense , did not associate with a diploid species, and 
therefore represented polyploid derivatives of diploid taxa that were not sampled and possibly 
are extinct. Associations of diploid and polyploid species as determined by rbcL analysis are 
consistent with morphological hypotheses of polyploid origins with the exception of B. ascendens, 
which did not possess the rbcL sequence of either of its putative diploid progenitors. 


Species of Botrychium subgenus Botrychium, commonly called moonworts, 
represent the smallest, most inconspicuous, and most species-rich subgroup 
of the genus. Although the mountains of western North America harbor the 
greatest number of species, subgenus Botrychium includes species that are 
distributed throughout most of the world (Clausen, 1938). However, the ranges 
of individual species vary remarkably. The most common species, Botrychium 
lunaria (L.) Sw., grows on nearly every continent, whereas B. gallicomontanum 
Farrar & Johnson-Groh is known only from a single locality. For many moon- 
wort species, current knowledge of their ranges may reflect more the extreme 
difficulty involved in finding the plants in the field rather than their true dis¬ 
tributions (Wagner and Wagner, 1994). 

A distinguishing feature of Ophioglossaceae, including subgenus Botry¬ 
chium, is that the entire aerial portion of the plant is composed of a single 
leaf divided into a fertile segment and sterile segment (Bower, 1926; Bierhorst, 
1971; Foster and Gifford, 1974). In general, a single leaf is produced per year, 
and it is typically small in size, ranging from 1-15 cm in height. The small 
size and simplicity of the leaf have limited the number of discrete characters 
available for the construction of classifications and made the systematic^ of 
subgenus Botrychium particularly problematic and controversial (Tryon and 
Tryon, 1982). 

As circumscribed by Clausen (1938), subgenus Botrychium contained 6 dis¬ 
tinct species (Table 1) and 15 varieties or subspecies. Over the last several 


376 


AMERICAN FERN JOURNAL: VOLUME 85 NUMBER 4 (1995) 


Table 1. Species of Botrychium subgenus Botrychium. Preliminary sectional assignments based on 
morphological and cytological studies of Wagner and Wagner (pers. comm.). The number after each 
species name indicates the documented ploidy level of the species (2x = diploid, 4x = tetraploid, 6x = 
hexaploid), with x = 45. Asterisks indicate species sampled in the present study. 


Clausen (1938) 

Wagner & Wagner 
(pers. comm.) 

Ploidy 

level 


Section Lanceolatum 


B. lanceolatum (J. Gmel.) Angstr. 

*B. lanceolatum 

2x 


Section Lunaria 


B. 1 unaria (L.) Sw. 

*B. lunaria 

2x 


*B. crenulatum W.H. Wagner 

2x 


B. pallidum W.H. Wagner 

2x 


*B. minganense Viet. 

4x 


4 B. spathulatum W.H. Wagner 

4x 


4 B. ascendens W.H. Wagner 

4x 


Section Matricariifolium 


B. matricariifolium A. Braun 

* B. matricariifolium 

4x 

B. boreale Milde 

B, boreale 

9 

• 


*B. acuminatum W.H. Wagner 

4x 


*B. hesperium (Maxon & Clausen) 

4x 


W.H. Wagner & Lellinger 



*B. echo W.H. Wagner 

4x 


*B. pedunculosum W.H. Wagner 

4x 


4 B. pinna turn St. John 

4x 


4 B. pseudopinnatum W.H. Wagner 

6x 


Section Simplex 


B. simplex E. Hitchc. 

*B. simplex 

2x 

B. pumicola Coville 

*B. pumicola 

2x 


*B. montanum W.H. Wagner 

2x 


B. mormo W.H. Wagner 

2x 


B. gallicomontanum Farrar & 

4x 


Johnson-Groh 



Sectional status undetermined 



4 B. campestre W.H. Wagner & Farrar 

2x 


4 B. line are W.H. Wagner 

? 

V 


*B. pa rad ox um W.H. Wagner 

4x 


4 B. Xwatertonense W.H. Wagner 

4x 


decades the number of species has increased to 24 as a result of the detailed 
morphological and cytological studies of W. H. and F. S. Wagner (Table 1). 
Wagner and Wagner and other workers have used increasingly subtle morpho¬ 
logical characters initially to recognize new species, many of which are poly¬ 
ploid (Wagner and Wagner, 1981, 1983b, 1986, 1990a, 1990b; Wagner et. al, 
1984; Farrar and Johnson-Groh, 1991). Through the use of comparative field 
techniques (Wagner and Wagner, 1983a), environmentally stable and presum¬ 
ably genetically fixed morphological characters have been identified and cor¬ 
related with cytological observations to yield current concepts of species de¬ 
limitations (Wagner, 1993). 
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One major consequence of the morphological simplicity characteristic of 
subgenus Botrychium is that many species are extraordinarily difficult to dis¬ 
tinguish. Because the morphological characters necessary to distinguish spe¬ 
cies are subtle, many species of subgenus Botrychium are characterized as 
“cryptic” species. Cryptic species were defined by Stebbins (1950) as, “. . . 
population systems which were believed to belong to the same species until 
genetic evidence showed the existence of isolating mechanisms separating 
them.” According to Paris et al. (1989), cryptic species typically: 1) are poorly 
differentiated morphologically, 2) represent distinct evolutionary lineages be¬ 
cause they are reproductively isolated, and 3) have historically been misinter¬ 
preted as members of a single species. Many species of subgenus Botrychium 
meet these criteria. First, they are extremely difficult to differentiate without 
considerable field experience. Small (<2cm) or immature individuals can be 
nearly impossible to identify for even the most experienced collectors. Second, 
the existence of sterile, putative hybrids with intermediate morphology and 
abortive spores has provided evidence of hybridization and reproductive iso¬ 
lation between species (Wagner, 1980, 1991; Wagner and Wagner, 1988; Wag¬ 
ner et al., 1985). Third, the majority of newly recognized species was separated 
from within a group previously recognized as a single species. For example, 
B. pirmatum H. St. John, B. pseudopinnatum W.H. Wagner, B. pedunculosum 
W.H. Wagner, B. hesperium (Maxon and R.T. Clausen) W.H. Wagner and Lel- 
linger, B. echo W.H. Wagner, and B. acuminatum W.H. Wagner were all con¬ 
sidered minor variants of B. matricariifolium A. Braun (Clausen, 1938), but 
are now recognized as species (Wagner and Wagner, 1981, 1983b, 1986, 1990b). 

Plant groups for which morphological distinctions among species are subtle 
are appropriate candidates for molecular analyses. Although plastid inheri¬ 
tance is untested in Ophioglossaceae, molecular characters derived from plas¬ 
tid DNA often have the advantage of uniparental inheritance, either paternally 
(gymnosperms: Neale et al., 1989; Neale and Sederoff, 1988, 1989) or mater¬ 
nally (angiosperms: Harris and Ingram, 1991; ferns: Gastony and Yatskievych, 
1992). Sequencing of the plastid gene rbcL has provided insight into higher 
order relationships among angiosperms (Chase et al., 1993), red algae (Fresh¬ 
water et al., 1994), and ferns (Hasebe et al., 1995). Because of the conservative 
nature of rbcL, most studies have focused on relationships at the rank of genus 
or higher. Although rbcL has proved useful to address species-level questions 
in a limited number of plant groups (Price and Palmer, 1993; Williams et al., 
1994), the extent to which rbcL can address species-level relationships remains 
largely unexplored. Furthermore, studies utilizing rbcL for phylogenetic re¬ 
construction typically rely on single collections to represent entire species, 
even when these species have broad geographic distributions, apparently as¬ 
suming that within-species polymorphism in rbcL is non-existent or that any 
polymorphism present will not influence phylogenetic inferences. 

Representative species of the subgenera Osmundopteris, Sceptridium, and 
Botrychium from genus Botrychium had a mean rbcL sequence divergence val¬ 
ue of 2.4% (Hauk, Parks, and Chase, unpubl. data), a value comparable tc that 
observed among some angiosperm genera (Les, 1994). Sequence divergence 
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between two species of subgenus Botrychium, B. lanceolcitum (}. Gmel.) Angstr. 
and B. simplex E. Hitchc., was 1.1%. These levels of divergence indicated that 
rbcL might be informative enough to address questions of phylogenetic rela¬ 
tionships among species or species groups within subgenus Botrychium. Ac¬ 
cordingly, three goals for the study were identified: 1) to investigate phyloge¬ 
netic relationships among diploid Botrychium species and species groups, 2) 
to evaluate relationships among diploid and polyploid taxa in an effort to 
identify diploid progenitors oi polyploid species, and 3) to assess within-spe- 
cies consistency of rbcL sequences across broad geographic distances. 


Materials and Methods 


DNA Extraction, Gene Amplification and Sequencing. —Plant material from 
20 of the 24 currently recognized species of subgenus Botrychium was col¬ 
lected and silica-gel dried. Material of four species, B. mormo W.H. Wagner, 
B. pallidum W.H. Wagner, B. boreale Milde, and B. gallicomontanum, was not 
available tor analysis. Table 2 contains collection sites and voucher informa¬ 
tion numbers. Total genomic DNA was extracted from approximately 0.03 to 
0.50 grams of silica gel dried tissue that was ground at 65°C in 2X CTAB buffer 
(2% hexadecyl trimethyl ammonium bromide in 100 mM Tris-HCl pH 8.0 with 
1.4 M NaGl and 20 mM EDTA) after Doyle and Doyle (1987), as modified by 
Wendel (1989) (20% polyvinylpyrrolidone, 5.0 mM ascorbic acid, 4.0 mM die- 
thyldithiocarbamic acid sodium salt). Treatment with chloroform: isoamyl al¬ 
cohol (24:1) removed nonpolar substances. Total DNA was purified by density 
gradient centrifugation in CsCl 2 -ethidium bromide. After dialysis, the DNA 
was precipitated (1.0 M NaOAc, 70% EtOH) according to Sambrook et al. 
(1989) and spooled out of solution with a glass hook, rinsed in 70% ethanol, 
air dried, and resuspended in TE buffer pH 8.0. Amplification of the rbcL gene 
by standard methods (Saiki et al., 1987) used a forward primer (5'-ATGTCAC- 

GAGAAAGAGAAACTAAAGGAAGT-3') that attached to the first 30 base pairs 
of the exon, and a reverse primer (5'-CTTCACAAGCAGCAGCTAGTTCAG- 
GACTCC-3') that began at the exon 3' position 1352. Amplified product s were 
purified using glass beads (Vogelstein and Gillespie, 1979). The purified prod¬ 
ucts were sequenced directly using the dideoxynucleotide method of Sanger 
(1981). Both strands were sequenced using a combination of eight internal 
sequencing primers designed specifically for Ophioglossaceae (primer se¬ 
quences available on request), plus the two amplification primers. A total of 
1321 base pairs was obtained for each taxon (1.351 minus the 30 base pairs of 
the forward primer). All substitutions unique to subgenus Botrychium species 
were recorded and verified by re-examination of the sequencing autoradio¬ 
graphs from both forward and reverse primers. 


Data Analysis.—A total of 26 complete sequences, representing 20 species of 
subgenus Botrychium, was analyzed. Pairwise percent sequence divergences 
were calculated by dividing the number oi nucleotide differences between two 
taxa by the total number of nucleotides sequenced (1321 for all sequences). 
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Table 2. Species of subgenus Botrychium used in rbcL analysis, collector and voucher citation, Gen- 
bank accession number and collection site. All Hauk vouchers are deposited at NCU. Wagner vojchers 
are deposited at MICH. 


Source/ Genbank 

Species voucher number Collection site 


Moran 5426 (MO) L40966 
Hauk 564 L40965 


Botrychium tun aria 
Botrychium lunaria 
Battych ium minganense 
Botrychium minganense 
Botrychium minganense 
Botrychium minganense 
Botrychium crenulatum 
Botrychium ascendens 

W' 

Botrychiuni spathulatum 
Botrych ium spath u latum 
Botrychium simplex 
Botrychium simplex 
Botrychium pumicola 
Botrych i um man tan um 

wr 

Botrychium lanceolatum 

m 

var. angustisegmentum 
Botrych ium matrica riifolium 
Botrychium hesperium 
Botrychium echo 
Botrychium acuminatum 
Botrychium pinnatum 
Botrych iun i pseudopinnaturn 
Botrychium pedunculosum 
Bot rye hi um pa rad ox um 
Botrychium Xwatertonense 
Botrychium lineare 
Botrych ium campestre 


Hauk 578 

L40970 

Hauk 566 

L40971 

Hauk 598 

L40968 

Hauk 584 

L40969 

Hauk 616 

L40959 

Hauk 529 

L40982 

Wagner 88036 

L40980 

Hauk 562 

L40979 

Hauk 619 

L40978 

Hauk 561 

L40977 

Hauk 618 

L40976 

Hauk 607 

L40916 

Hauk 571 

L40963 

Hauk 551 

L40967 

Hauk 552 

L40960 

Hauk 595 

L40962 

Hauk 553 

L40922 

Hauk 604 

L40974 

Wagner 88037 

L40975 

Hauk 615 

L40973 

Hauk 610 

L40972 

Hauk 611 

L40981 

Hauk 581 

L40964 

Farrar s.n. (ISU) 

1.40961 


Taiwan 

Marathon, Ontario 
Grand Sable Dunes, MI 
Pic River, Ontario 
Boulder Co., CO 
Independence Pass, CO 
Hurricane Creek, OR 
Hurricane Creek, OR 
Angler Settlement, Ontario 
Marathon, Ontario 
Mt. Ashland, OR 
Grand Sable Dunes, Ml 
Newberry Caldera, OR 
Lake Co., MT 
Taquamenon Falls, MI 

Grand Sable Dunes, MI 
Grand Sable Dunes, MI 
Echo Lake, CO 
Grand Sable Dunes, MI 
Stagger Inn. WA 
Angler Settlement, Ontario 
Lostine River, OR 
Waterton, Alberta 
Waterton, Alberta 
Pike’s Peak, CO 
Beamis Creek, IA 


For cladistic analyses, B. multifidum (J. Gmel.) Rupr. and B. lunaricides 
(Michx.) Sw. of subgenus Sceptridium were chosen as outgroups, based on a 
larger rbcL analysis of Ophioglossaceae (Hauk, Parks, and Chase, unpubl. data). 
All sequences were analyzed using PAUP (Phylogenetic Analysis Using Par¬ 
simony) version 3.1.1 (Swofford and Begle, 1993). Fitch (1971; equal weights) 
parsimony analysis of 1000 random taxon entries was performed using Tree 
Bisection-Reconnection (TBR) swapping with MULPARS selected. One anal¬ 
ysis was performed on rbcL sequences of diploid species and a second analysis 
was conducted with sequences of both diploid and polyploid species. The 
chromosome number for B. lineare W.H. Wagner is not known, but for pur¬ 
poses of analysis it was assumed to be diploid, based on isozyme banding 
patterns (Hauk, 1990). Bremer (1988) support (decay of parsimony) identified 
the tree lengths at which each branch collapsed into a polytomy. Figure 1 
reports Bremer support values as a “d” preceeding a number (i.e., “d4” means 
that the branch decayed at four steps shorter than the length of the single most 
parsimonious tree). As a comparison to heuristic methods of analysis, I con- 
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^ 1 

diploid species of Botrychium subgenus Botrychium. All characters were equally weighted and 
unordered. Sequences ol B. multifidum and B. lunarioides (subgenus Sceptridium) were used as 
outgroups. Segment lengths are listed above the branches, and the number of steps of relaxed 
parsimony required to collapse a branch are listed below the branches. The tree has a length of 
40 steps, 0 = 0.918, RI=0.930, and HI = 0.083. 


ducted an exhaustive search of the 9 different rbcL sequences represented 
among the 20 species (i.e., some species had identical sequences). 


Results and Discussion 


Levels of Sequence Variation in rbcL. —Analysis of all sequences identified 
six species groups with identical rbcL sequences: 1) B. lanceolatum (var. an- 


Moore) 


ifoli 


B. 


W.H. W 


W.H. Wagner, B. Xwntertonense W.H. W 


illations of B. mingancnse Viet.; 4) B. simplex (MI) and B. montanum W 
Wagner; 5) B. simplex (OR) and B. pumicola Cov.; and 6) B. campestre W 
Wagner and Farrar and B. spathulatum W.H. Wagner. Botrychium lunnria fil 
wan) and B. lineare W.H. Waener had the onlv two rbcf, senuences that m 
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sessecl autapomorphic mutations. The mean ot pairwise comparisons ot se¬ 
quence divergence among the 9 different rbcL sequences distributed among 20 
subgenus Botrychium species was 0.6%, a value much lower than the 1.87% 
reported for sister species of Polypodium (Haufler and Ranker, 1995). However, 
diploids representing each of the Wagner and Wagner sections had a mean 
sequence divergence of 1.1%, a level comparable to that observed among gen¬ 
era of aquatic cresses (Les, 1994). Levels of sequence divergence at a given 
rank will undoubtedly differ greatly among different groups of plants, due to 
different rates of mutation or contrasting concepts of taxonomy. Although the 
0.6% mean sequence divergence among subgenus Botrychium species repre¬ 
sents a very low level of divergence, statistics such as percentage sequence 
divergence do not in themselves provide a reliable evaluation of the infor¬ 
mation content of the divergence observed, or an indication of the strength of 
analyses based on them. Ritland and Eckenwalder (1992) estimated that for 
phylogenetic analyses optimal levels of sequence divergence occur around 
10%, a point at which the number of single-hit substitution sites is maximized 
relative to the number of multiple-hit substitution sites. High levels of se¬ 
quence divergence may lead to relationships that have poor internal and ex¬ 
ternal support (Manhart, 1994). At lower levels of sequence divergence 
(< 10%), the principle constraint for phylogenetic studies is identifying 
enough synapomorphies to resolve relationships (Les, 1994). 

The absence of sequence divergence among many species of subgenus Bo¬ 
trychium indicates that rbcL is too conservative to address all species-level 


variation detected among 


informative 


homoplasy (see below). The conservative nature of rbcL coupled with its pu¬ 
tative uniparental pattern of inheritance should reveal clear patterns of rela¬ 
tionships, first, among diploid species groups and, second, between polyploids 
and their diploid chloroplast progenitor. 

Phylogenetic Analyses. —Forty-three variable nucleotide sites were detected 
among the 22 species included in this study. Of these 43 variable sites, 24 
were variable within subeenus Botrychium. Sixteen sites were informative in 


informative 


parsimony 


most parsimonious 


consistency index (Cl) of 0.918, a retention index (RI) of 0.930 and a homo¬ 
plasy index (HI) of 0.082 when both informative and uninformative characters 
were used. When only phylogenetically informative characters were included, 
a tree of 34 steps resulted with a Cl of 0.882, a RI of 0.930 and a HI of 0.118. 
In the second analysis, which included sequences of both diploid and poly¬ 
ploid species, a single most parsimonious tree of 49 steps was produced (Fig. 
2) with a Cl of 0.918, a RI of 0.974 and a HI of 0.082. When uninformative 
characters were excluded the tree had 41 steps, a Cl of 0.902, a RI of 0.974 


search vielded the same 


taxa 
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G 


0 


0 


0 


0 


1 


0 


0 


0 


1 


0 


0 


G 


G 


0 


0 


G 


G 


G 


G 


1 


£ 


B. spathulatum 

4x 

Ontario 

B. spathulatum 

4x 

Ontario 

B. lineare 

?x 

Colorado 

B. campestre 

2x 

Iowa 

B. ascendens 

4x 

Oregon 

B. simplex 

2x 

Oregon 

B. pumicola 

2x 

Oregon 

B. simplex 

2x 

Michigan 

B. montanum 

2x 

Montana 

— B. paradoxum 

4x 

Alberta 

B. Xwatertonense 

4x 

Alberta 

B. minganense 

4x 

Michigan 

B. minganense 

4x 

Colorado 

B. minganense 

4x 

Colorado 

— B, minganense 

4x 

Ontario 

B. lunaria 

2x 

Taiwan 

B. lunaria 

2x 

Ontario 

B. crenulatum 

2x 

Oregon 

B. lanceolatum 

2x 

Michigan 

B. hesperium 

4x 

Michigan 

B. matricariifolium 

4x 

Michigan 

B. pseudopinnatum 

6x 

Ontario 

B.echo 

4x 

Colorado 

B. pedunculosum 

4x 

Oregon 

B. pinnatum 

4x 

Washington 

B. acuminatum 

4x 

Michigan 

B. multifidum 



B. lunarioides 




Fig. 2. The single most parsimonious Fitch tree from the analysis of 26 rbcL sequences from 20 
species of Botry'chium subgenus Botrychium. All characters were equally weighted and unordered. 
Sequences of B. multifidum and B. lunarioides (subgenus Sceptridium) were used as outgroups. 
Segment lengths are listed above the branches, and the number of steps of relaxed parsimony 
required to collapse a branch are listed below the branches. The tree has a length of 49 steps, 
CI = 0.918, Rl = 0.974, and HI = 0.082. • 


Within subgenus Botrychium there are three well-supported diploid species 
clades (Fig. 1): The lanceolatum clade (one species), the lunaria clade (two 
species), and the simplex-campestre clade (five species). The lunaria and sim- 
plex-campestre clades collapse on a consensus tree four or more steps longer 
than the shortest tree (d4). Within the simplex-campestre clade there were two 
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less well-supported subclades: The simplex subclade (three species, dl) and 
the campestre subclade (two species, d2). Within the simplex subclade, two 
populations of B. simplex differed by a single substitution. Botrychium sim¬ 
plex from Oregon and B. pumicola had identical sequences, as did B. simplex 
from Michigan and B. montanum from Montana. Two species, B. lunaria and 
B. crenulatum, formed the lunaria clade. 

The second analysis included sequences of both diploid and polyploid 
subgenus Botrychium species (Fig. 2). Most rbcL sequences of polyploid 
species were identical to that of one of the diploid species. The polyploid 
species, B. pinnatum, B. pseudopinnatum, B. matricariifolium, B. hesper- 
ium, B. echo, B. pedunculosum, and B. acuminatum, had rbcL sequences 
identical to the diploid B. lanceolatum, and both populations of the tetra- 
ploid species, B. spathulatum, shared the same rbcL sequence with B. cam¬ 
pestre. Four polyploid subgenus Botrychium species did not share identical 
rbcL sequences with any diploid species sampled. The tetraploid species, 
B. ascendens W.H. Wagner, associated with the campestre subclade, but 
lacked one of the two synapomorphies that defined the subclade. Three 
other taxa, B. minganense, B. paradoxum, and B. Xwatertonense did not 
possess synapomorphies of the campestre or simplex subclades, but nev¬ 
ertheless possessed the four synapomorphies that defined the simplex-cam- 
pestre clade. The analysis detected no autapomorphic substitutions that dis¬ 
tinguished B. paradoxum from B. Xwatertonense or among the four popu¬ 
lations of B. minganense. 


The Monophyly of Subgenus Botrychium. —Species of subgenus Botry¬ 
chium possess a suite of morphological and anatomical characters that dis¬ 
tinguish them collectively from other subgenera of genus Botrychium (Bow¬ 
er, 1926; Clausen, 1938; Tryon and Tryon, 1982). The small size, pale color, 
fleshy leaves, pinnately compound trophophores (except B. lanceolatum, 
see below), invariable presence of a fertile segment, and ephemeral leaf 
growth support the monophyly of the subgenus. Micromorphological char¬ 
acters, such as similar pollen architecture (Tryon and Lugardon, 1991) and 
a leaf sheath covering with a fused slit (Kato, 1987), have been reported as 
synapomorphies for the group. Reports of hybrids among many species doc¬ 
ument the close genetic relationship among species (Wagner, 1980, 1991; 
Wagner et al., 1985; Wagner and Wagner, 1988). No reports of hybrids be¬ 
tween members of subgenus Botrychium and the often sympatric subgenera 
Sceptridium and Osmundopteris are known, and further suggest isolating 
mechanisms among the subgenera. 

In an rbcL analysis of Ophioglossaceae, all species ol subgenus Botry¬ 
chium form a monophyletic group (Hauk, Parks and Chase, unpubl. data). 
The addition of various outgroups within and outside of Ophioglossaceae 
does not change the strong support for the monophyly of subgenus Botry¬ 
chium. Although no morphological data are included in the analysis, all 
evidence, morphological and molecular, is consistent with the hypothesis 
that subgenus Botrychium is monophyletic. 
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Relationships Among Diploid Species of Subgenijs Botrychium. —Of the 

three major clades in the rbcL analysis, the lanceolatum clade is sister to 
the rest of the subgenus. 'The lunaria cdade (two diploid species) is sister to 
the simplex-campestre clade (the remaining five diploid species). Relation¬ 
ships among these three major clades are consistent with ideas of relation¬ 
ship based on morphology. Diploid B. lanceolatum is the only species in 
subgenus Botrychium that exhibits ternate division of the sterile segment, 
as do species of subgenera Sceptridium, Japanobotrychium and Osmundop- 
teris. If the presence of ternate leaf organization in B. lanceolatum is an¬ 
cestral and the pinnate leaf organization of the remaining species is derived, 
then morphological patterns are consistent with the rbcL data. 

The lunaria clade contains the diploid species B. lunaria and B. crenulatum, 
which are morphologically similar due to oppositely arranged, broadly fan¬ 
shaped pinnae. The lunaria clade is sister to the simplex-campestre clade. 
Within the latter, subclade simplex contains the diploids B. simplex, B. pum- 
icola and B. montanum, all of which share a cup-like curve in the terminal 
pinna of the sterile segment. The campestre subclade contains the diploid spe¬ 
cies B. campestre and the recently described B. lineare, for which ploidv is 
not known. A close relationship between B. campestre and B. lineare is evi¬ 
dent from morphology (Wagner, 1994), and the two species differ by only a 
single autapomorphic substitution in B. lineare. Thus, each of the diploid cla¬ 
des identified by rbcL analysis contains species that share similar morpholog¬ 
ical features. However, fundamental phylogenetic relationships among these 
diploid clades are difficult to establish based on morphological characters, 
because few characters exist and, of those, few are informative regarding re¬ 
lationships among clades. 

Studies of isozyme variation among closely related species are useful to es¬ 
tablish basic ideas of relationship (Crawford, 1983; Werth, 1989), although 
these data usually have been analyzed phenetically, and therefore do not spe¬ 
cifically address cladistic relationships. Nonetheless, hypotheses of relation¬ 
ships based on isozyme genetic identities (Nei, 1972) of six diploid species of 
subgenus Botrychium (Hauk, 1990) corroborate ideas of relationships based on 
parsimony analysis of rbcL sequences (Fig. 3). Diploid species with similar 
rbcL sequences tend to have higher isozymic genetic identities than diploids 
with divergent rbcL sequences. For example, diploid species B. crenulatum 
shares a genetic identity ol 0.53 with B. lunaria, the only other member of its 
rbcL clade. Species from other rbcL clades, such as B. simplex (0.34) and B. 
lanceolatum (0.39), have lower genetic identities with B. crenulatum than does 
B. lunaria. Botrychium lineare, a member of the simplex-campestre clade, has 
higher identities with B. simplex (0.59) and B. pumicola (0.83) than it does to 
species in other rbcL clades such as B. lunaria (0.24) and B. lanceolatum (0.04). 
The degree of correlation between isozyme and rbcL data supports the asser¬ 
tion that rbcL is informative concerning phylogenetic relationships of species 


groups in subgenus Botrychium. 


Intraspecific Variation in rbcL Sequences of Diploid Species. —Because most 
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B. lunaria 


B. crenulatum 


B. lineare 


-B. pumicola 


B. simplex 


B. lanceolatum 



Fic;. 3. Cluster analysis of six diploid species of Botrychium subgenus Botrychium. based on Nei's 
(1972) genetic identity of isozyme data, using the unweighted pair-group method (IJPGMA). The 
scale represents the range of potential genetic identity values. 


studies utilizing rbcL sequences focus on higher order relationships among 
genera and families, sometimes one accession of a single individual is assumed 
to represent a species (or even a genus or a family). At low levels of relation¬ 
ship (i.e., among species), where sequence divergence is minimal, within-spe- 
cies variation in rbcL may confuse phylogenetic relationships. One conse¬ 
quence of within-species rbcL sequence variation may be fixation of different 
base pair substitutions in different populations across the range of the species. 
Although few studies have investigated the possibility of within-species vari¬ 
ation in rbcL, species with broad distributions have the potential to possess 
more than one rbcL sequence. 

Two geographically widespread diploid species from subgenus Botrychium 
were selected for intraspecific sampling: B. simplex and B. lunaria. Intraspe¬ 
cific morphological variation is evident in B. simplex in which individuals 
collected at a single site may show considerable variation in the size and com¬ 
plexity of the lowermost pinnae pair. Some B. simplex individuals have large, 
well developed basal pinnae that are as dissected as the main axis of the sterile 
segment, whereas other individuals show no lobing or enlargement of lower 
pinnae. Historically, eastern and western varieties of B. simplex have been 
recognized (Clausen, 1938; Wagner and Wagner, 1993). For this study one pop¬ 
ulation from Michigan represents the eastern form and a population from Or¬ 
egon represents the western form. Botrychium lunaria is the most geographi¬ 
cally widespread species in subgenus Botrychium; collections are reported 
from nearly every continent (Clausen, 1938). The potential for long distance 
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dispersal is high in this species. Two widely disjunct populations, one from 
Michigan and one from Taiwan, represent the range of B. lunaria. 

Botrychium lunaria from Ontario does not possess the autapomorphic sub¬ 
stitution that distinguishes the Taiwan population of B. lunaria from B. cren- 
ulatum. If within-species sampling had not been conducted, then B. lunaria 
and B. crenulatuin would appear differentiated by a single substitution. Thus, 
the North American collections of B. lunaria and B. crenulatuin are more sim¬ 
ilar to each other than the two populations of B. lunaria collected from widely 
separate geographic areas, Ontario and Taiwan. Several explanations for this 
pattern of relationship are possible: 1) B. crenulatuin and B. lunaria may share 
a common ancestor that gave rise to geographically proximal, but morpholog¬ 
ically distinct species with identical rbcL sequences, whereas geographically 
disjunct populations of B. lunaria are morphologically indistinguishable but 
have divergent rbcL sequences: 2) introgression of the lunaria chloroplast type 
into B. crenulatuin may explain the lack of rbcL sequence divergence between 
B. lunaria and B. crenulatunr, 3) lineage sorting of chloroplast DNA polymor¬ 
phisms could produce the observed pattern of relationship between B. lunaria 
and B. crenulatuin (Doyle et ah, 1990). According to the lineage sorting hy¬ 
pothesis, rbcL polymorphisms were present in the common ancestor of B. lu¬ 
naria and B. crenulatuin, and both species inherited the polymorphisms. The 
Ontario population of B. lunaria became fixed for one rbcL pattern and the 
Tiawan population for a second rbcL pattern. By chance, B. crenulatuin be¬ 
came fixed for the same rbcL pattern as the Ontario population of B. lunaria, 
thus giving the appearance of most recent shared ancestry between B. crenu¬ 
latuin and the Michigan population of B. lunaria. 

Relationships among B. pumicola, B. montanum, and the two populations 
ol B. simplex are more complex. The western population of B. simplex shares 
a substitution with the western species, B. pumicola, whereas the eastern pop¬ 
ulation ol B. simplex lacks this single substitution, as does B. montanum, an¬ 
other western North American species. The rbcL tree indicates that the com¬ 
mon ancestor of the simplex subclade gave rise to B. montanum and B. simplex 
directly. Eastern and western populations of B. simplex then diverged, and the 
western population produced B. pumicola. Alternatively, phylogenetic sorting 
(Doyle et ah, 1990) or introgression (Rieseberg, et ah, 1990; Rieseberg and 
Brunsfeld, 1991) may account for the same patterns of relationship, although 
without further information from nuclear markers this hypothesis cannot be 
evaluated. 

The intraspecific variation detected in B. lunaria and B. simplex signifi¬ 
cantly influences the interpretation of the relationship of these two species to 
that ol closely related species. In both B. simplex and B. lunaria, sampling of 
a single population would have resulted in erroneous interpretations of species 
relationships. Furthermore, these data suggest the possibility of undetected 
rbcL variation in these and other diploid subgenus Botrychium species. 


Associations of Diploid and Polyploid Species. —The ancestry of many pu¬ 
tative allopolyploid species of subgenus Botrychium is enigmatic, because 
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Table 3. Hypothetical parents for polyploid species of subgenus Botrychium (after F.S. Wagner, 1993). 


Polyploid species 


Hypothetical parents 


B. minganense 
B. spathulatum 
B. a seen dens 
B. gallicomontanum 
B. pinnatwn 
B. hesperium 
B. echo 

B. pedunculosum 
B. matricarii folium 
B. acuminatum 
B. pseudopinna turn 

B. Xwatertonense 


B. lunaria X B. pallidum 
B. lunaria X B. campestre 
B. crenulatum X B. montanum 


B. simplex X B. campestre 
B. lanceolatum X B. lunaria 

B. lanceolatum X B. simplex 

B. lanceolatum X B. campestre 

B . lanceolatum X B. montanum 

B . lanc eolatum X B. pallidum 

metaspecies of B. matricariifoli 
B. pinnatum X B. simplex 


urn 




there are relatively few discrete morphological characters that distinguish taxa. 
Hybridization followed by polyploidy blurs species specific characters and 
makes the recognition and identification of species difficult. Assuming uni¬ 
parental inheritance of the chloroplast genome, recent polyploid derivatives 
of extant diploid species should possess similar or identical rbcL sequences to 
known diploids. Polyploid species of subgenus Botrychium associate with var¬ 
ious diploid species in the rbcL tree (Fig. 2) and these associations are largely 
congruent with previously hypothesized ideas of polyploid origin (Wagner, 
1993; Table 3). 

The lanceolatum clade contains seven polyploid species and a single dip¬ 
loid, B. lanceolatum. The seven polyploid species have rbcL sequences iden¬ 
tical to that of B. lanceolatum and, based on morphological studies, are pu¬ 
tative allopolyploid derivatives of B. lanceolatum and pinnate diploids (Wag¬ 
ner, 1993). Thus, B. lanceolatum appears to be the chloroplast parent of the 
seven polyploid species in the lanceolatum clade. Sequence divergence in rbcL 
almost certainly does not coincide with all speciation events, and is thus not 
useful to distinguish all subgenus Botrychium species. A relatively recent or¬ 
igin for each of the polyploid species in the lanceolatum clade could explain 
why none of the seven polyploids have autapomorphic substitutions and yet 
all possess the six substitutions that distinguish B. lanceolatum from other 
diploid species. 

The second major rbcL clade is the lunaria clade, which includes only the 
diploids B. lunaria and B. crenulatum. Botrychium crenulatum was the only 
diploid in addition to B. lunaria placed in section Lunaria by Wagner and 
Wagner (pers. comm.). Although morphological evidence exists for B. lunaria 
as one diploid progenitor of B. spathulatum, B. minganense, B. pinnatum , and 
B. pseudopinnatum (via B. pinnatum), none of the polyploids possess the rbcL 
sequence of B. lunaria. Because B. lunaria is the most geographically wide¬ 
spread of all subgenus Botrychium species, it is surprising that no polyploids 
share its rbcL sequence, but further sampling among putative polyploid deriv- 
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atives of B. lunaria might uncover examples of inheritance of the lunaria rbcL 
sequence. 

The simplex-campestre clade contains ten species, five of which are poly¬ 
ploids. The clade is composed of three distinct groups, two poorly supported 
subclades (dl) and three taxa that did not associate with either subclade. The 
campestre subclade contains the diploid B. campestre, the tetraploids B. spa- 
thulatum and B. ascendens, and one additional species whose ploidy level is 
not known, B. lineare. All four species exhibit the common feature of narrow, 
almost linear pinnae, with the possible exception of B. spathulatum, in which 
the broader pinnae shape is probably due to the influence of B. lunaria as one 
of its diploid parents (Wagner, 1993). The identical rbcL sequences of B. spa¬ 
thulatum and B. campestre implicate B. campestre as the chloroplast parent 
oi B. spathulatum. In this case the more geographically restricted and less 

common diploid parent, B. campestre, contributed its chloroplast genome to 
the polyploid. 

The association ol B. ascendens with the campestre subclade was not sug¬ 
gested in previous morphological investigations (Wagner, 1993; Table 3), al¬ 
though all members ol the subclade possess narrow pinnae. Botrychium as¬ 
cendens does not possess both of the synapomorphies of the campestre sub¬ 
clade. The lack of one synapomorphy may be due to early divergence of B. 
ascendens from the common ancestor of B. campestre, B. lineare, and B. spa¬ 
thulatum. Alternatively, an unsampled diploid with the B. ascendens chloro¬ 
plast type may exist. 

The tetraploids, B. mmganense, B. paradoxum and B. Xwatertonense do not 
possess the synapomorphies that distinguish the simplex and campestre sub¬ 
clades and appear to retain the ancestral rbcL sequence of the entire simplex- 
campestre clade. There are at least three explanations for the lack of associa¬ 
tion of these tetraploid species with any known diploid: 1) B. minganense, B. 
paradoxum and B. Xwatertonense have reversions of synapomorphies that es¬ 
tablish the campestre and simplex subclades as distinct (an explanation that 
seems unlikely given the conservative nature of rbcL)\ 2) B. minganense, B. 
paradoxum, and B. Xwatertonense may have diverged from the common an¬ 
cestor of the simplex-campestre clade, and thus do not possess the synapo¬ 
morphies of either the simplex or campestre subclades; 3) B. minganense, B. 
paradoxum and B. Xwatertonense may possess the chloroplast genome of a 
diploid that was not sampled in this analysis [B. pallidum) or is extinct. Ma¬ 


terial ot the diploid B. pallidum, a potential member of this clade, should be 
available in the near future. Only extensive collections and continued field¬ 
work by experienced collectors can address the possibility of the existence of 
an unknown diploid. 

The patterns of relationship between diploids and polyploids reveal several 
interesting speciation trends in subgenus Botrychium. First, certain diploids 
[B. lanceolatum and B. campestre) are frequently involved as the chloroplast 
parent in hybridization events leading to polyploidy. Second, species such as 
B. lunaria, B. crenulatum, and three species of the simplex subclade ( B. sim¬ 
plex, B. montanum, B. pumicola ) do not appear to be chloroplast parents of 
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any polyploids sampled. Members of the lunaria clade and the simplex sub- 
clade appear to be speciating at the diploid level while presumably contrib¬ 
uting only nuclear genomes to the polyploids. Third, the campestre subclade 
is the only group containing both closely related diploid species and a poly¬ 
ploid derivative possessing the rbcL pattern of a diploid subclade member. 
Reasons for these patterns of relationship are enigmatic. Differential produc¬ 
tion of antheridia and archegonia in gametophytes of different diploid species 
could account for these patterns, although no data concerning overall produc¬ 
tion of antheridia and archegonia in gametophytes of subgenus Botrychium are 
available at this time. 


Intra,specific Variation in Polyploid Species. —Within-species plastid DNA 
variation in allopolyploid taxa may originate from reciprocal crosses when 
there are multiple hybridization events between progenitor diploids (Gibby, 
1977; Werth et ah, 1985; Haufler and Soltis, 1986; Barrington et al., 1989; Soltis 
and Soltis, 1989; Stein and Barrington, 1990). Different parents may contribute 
their chloroplast genome to the resulting hybrids. Thus, conspecific popula¬ 
tions of allopolyploids may possess both parental chloroplast genomes (Stein 
and Barrington, 1990). In the tetraploid, B. minganertse, a high level of mor¬ 
phological variation may indicate the possibility of unrecognized cryptic spe¬ 
cies (W. H. Wagner, pers. comm.). High amounts of isozymic variation among 
eleven B. mirtganense populations lends support to the hypothesis of a com¬ 
plex origin for B. minganense (Hauk, 1990). In order to test for different chlo¬ 
roplast genomes in B. minganense, four populations, two from the Great Lakes 


sam 


re- 


lsozymic 


four populations of B. minganense share identical rbcL sequences. This lack 
of rbcL variation is consistent with the hypothesis of a single plastid origin for 
eastern and western B. minganense populations. 

Bimodal morphological variation is apparent in the putative allotetraploid, 
B. spathulatum ; different individuals may strongly resemble one or the other 
of the putative parents, B. campestre or B. lunaria (Hauk, unpubl. data), per¬ 
haps as a result of introgression. Two populations of B. spathulatum were 
collected from southern Ontario, one which more closely resembled B. cam¬ 
pestre and the other which strongly resembled B. lunaria. Sequences from both 
B. spathulatum populations are identical to the rbcL sequence of B. campestre, 
and are consistent with a single origin for these populations of B. spathulatum. 
The lack of rbcL variation in B. minganense and B. spathulatum indicates that 
hybridization events leading to allopolyploid speciation are extremely rare 
and/or consistently favor one parent as the chloroplast donor. 


Comparisons of bbcL Tree and Morphological Ideas of Relationship. —The 

Wagner and Wagner (pers. comm.) classification of subgenus Botrychium dif¬ 
fers from the rbcL tree in four ways (Table 4): 1) the rbcL tree fuses the Wagner 
and Wagner sections Lanceolatum and Matricariifolium ; 2) the rbcL tree sup¬ 
ports the establishment of a new section of subgenus Botrychium, provision¬ 
ally referred to as “section Campestre” pending further morphological study 
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Tablk 4. Contrasting sectional classifications for species of Botrychium subgenus Botrychium , based 
on Wagner and Wagner's (pers. comm.) most recent, as yet not fully published system and the present 
rbcL analysis. Species not sampled for the rbcL analysis are indicated with a dash (—). Note that “section 
Canipestre" is a provisional name. 


Species 

Wagner and Wagner 

C L? 

(sectional placement) 

rbcL analysis 

«r 

(sectional placement) 

B. lanceolatum 

Lanceolatum 

Lanceolatum 

B. matricariifolium 

Matricariifolium 

Lanceolatum 

B. acuminatum 

Matricariifolium 

Lanceolatum 

B. hesperium 

Matricariifolium 

Lanceolatum 

B. echo 

Matricariifolium 

Lanceolatum 

B. pedunculosum 

Matricariifolium 

Lanceolatum 

B. pinnatum 

Matricariifolium 

Lanceolatum 

B. pseudopinnatum 

Matricariifolium 

Lanceolatum 

B . boreale 

Matricariifolium 

4/ 


B. simplex 

Simplex 

Simplex 

B. pumicola 

Simplex 

Simplex 

B. montanum 

Simplex 

Simplex 

B. mormo 

Simplex 


B. gallicomontanum 

Simplex 


B. lnnaria 

Lnnaria 

Lnnaria 

B. crennlatnm 

Lnnaria 

Lnnaria 

B. pallidum 

Lnnaria 


B. ascendens 

Lnnaria 

“ Canipestre ” 

B. spathulatum 

Lnnaria 

“ Canipestre " 

B. minganense 

Lnnaria 

Status undetermined 

B. paradoxum 

Status undetermined 

Status undetermined 

B. Xwatertonense 

Status undetermined 

Status undetermined 

B. lineare 

Status undetermined 

“ Canipestre" 

B. canipestre 

Status undetermined 

“ Canipestre " 


and characterization, which would include B. canipestre, B. lineare, and B. 
spathulatunr, 3) B. ascendens aligns wilh “section Canipestre" and not with 
one of its hypothetical diploid parents, B. crennlatnm or B. montanum ; 4) B. 
minganense does not associate with members of section Lnnaria in the rbcL 
tree, and its sectional affinity is unknown, although it groups closely with B. 
paradoxum and B. Xwatertonense. 

Despite these differences, the rbcL tree of subgenus Botrychium reflects to a 
degree the relationships as hypothesized by the morphological studies of Wag¬ 
ner and Wagner. The major discrepancies between the two classifications may 
be attributed to problems associated with interpretation of biparentally inher¬ 
ited morphological characters versus uniparentallv inherited rbcL sequences. 
For example, the rbcL tree places B. ascendens, B. spathulatum, and B. min- 
ganense outside section Lnnaria. although either B. lnnaria or B. crennlatnm 
is probably involved in the parentage of these polyploid species, but not as 
the chloroplast donor. Morphological phenotypes of allopolyploids would not 
be expected to more closely resemble the chloroplast parent than the non- 
chloroplast parent, i.e., the allopolyploids would be placed closest to the par- 
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ent to which the closest morphological resemblance occurs, whether or not it 
is the chloroplast donor. 

Classification of hybrids and allopolyploid taxa is particularly problematic 
if the parental taxa belong to different taxonomic groups (Kellogg, 1989; Wag¬ 
ner, 1983), because an allopolyploid species does not belong exclusively to 
the taxonomic group of either parent. Debate concerning the treatment of hy¬ 
brids and allopolyploids in classifications has produced at least three possible 
solutions to the dilemma (Wagner, 1983): a) placement of hybrid or allopoly¬ 
ploid species next to the first parent appearing in the classification and ref¬ 
erence to the hybrid or allopolyploid again when the second parent is listed; 
b) placement of all hybrids and allopolyploids in a special section at the end 
of the species treatment; or c) arbitrary placement of hybrids and allopoly¬ 
ploids either with the first-listed parent or in a special section, based on char¬ 
acters such as fertility versus sterility, low versus high abundance, or other 
factors deemed justifiable by the researcher. 

These rbcL data provide strong evidence of the probable chloroplast parent 
of many subgenus Botrychium polyploids, although identification of the non- 
chloroplast parent(s) remains uncertain for many species, due to the inherent 
limitations of interpreting cryptic morphological differences that are further 
obscured by hybridization events. Other subgenus Botrychium species {B. nun- 
garxense, B. paradoxum, and B. ascendens ) do not have demonstrable chlo¬ 
roplast progenitors, and their origin remains a mystery. The most appropriate 
solution to the problem of classifying subgenus Botrychium polyploids is the 
creation of a special section in which all putative allopolyploids are placed. 
The allopolyploid section should be retained until more rigorous evidence 
elucidates the non-chloroplast progenitor of these species. Only when a clear 
understanding of the ancestry of these species is available will it be possible 
to create a classification system that reflects the evolutionary history of allo¬ 
polyploid species of subgenus Botrychium. 


Conclusions 

Even low amounts of sequence divergence among species (0.6% in subgenus 

Botrychium ) can contain phylogenetically informative signal if coupled with 

low levels of homoplasy. The divergence detected among species of subgenus 

Botrychium demonstrates that rbcL is more applicable to species-level ques- 
— 

tions than previously thought, although certainly not all species-level distinc¬ 
tions may be resolved. In the absence of more variable yet easily studied genes, 
rbcL can provide basic evidence concerning phylogenetic relationships among 
sneeies and snecies erouDS. Within submenus Botrychium, analysis of rbcL data 


provides the first opportunity to establish fundamental phylogenetic relation¬ 
ships among a confusing assemblage of morphologically cryptic species. In 
addition to supporting the monophyly of subgenus Botrychium, the rbcL phy- 
logeny identified B. Icmceolatum as sister to all other diploid species, and the 
group composed of B. lunaria and B. crenulatum as sister to the five diploid 
species of the simplex-campestre clade {B. campestre, B. linecire, B. simplex. 
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B. montanum, and B. pumicola). The rbcL phylogeny supports the establish¬ 
ment of a new section (provisionally called “section Campestre”) containing 
B. campestre, B. lineare, and B. spathulatum. Relationships based on rbcL 
analyses are largely consistent with morphological and isozymic hypotheses 
of relationship. 

In groups possessing polyploids of uncertain origin, rbcL can provide strong 
evidence of the chloroplast parent of each polyploid, provided that the dip¬ 
loids are well differentiated. Analysis of rbcL sequences establishes affinities 
of polyploid species of subgenus Botrychium to their diploid chloroplast do¬ 
nors, and thus provides the strongest evidence to date concerning polyploid 
origins. The within-species rbcL variation detected in B. lunaria and B. sim¬ 
plex demonstrates that intraspecific rbcL variation can significantly influence 
interpretations of species relationships. Population sampling for studies of spe¬ 
cies relationships may be necessary to avoid misinterpretation of phylogenetic 
relationships. 
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